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I. INTRODUCTION. 


Tue ordinary theory of errors —the method of least squares — has 
been developed in three ways (methods of Gauss, Hagen, and Crofton) ; 
not only is each of these developments based upon assumptions to which 
not every one will assent, but at certain stages in each of the two later 
developments exactness is replaced by approximation through the neglect 
of certain quantities (or terms) regarded as very small (or productive of 
very little effect) in comparison with others retained, whereas this neglect 
is not fully justified, as it seems to us. 

The present paper is the outcome of an attempt to construct a theory 
of errors upon such simple principles as will be generally admitted to be 
necessary for the mathematical treatment of any such theory. The im- 
mediate object of the theory of errors is to determine the probability that 
the error made in the direct measurement of a physical magnitude shall 
lie between given limits. More complex applications of the theory 
(e. g- to functions of quantities directly observed and to quantities deter- 
mined by calculations based upon direct observations) can be connected 
with the solution of the fundamental problem just stated by means 
familiar to all who have used the ordinary theory. Of course we con- 
sider only accidental errors, assuming the observations to have been 
corrected for all systematic errors before being subjected to our investi- 
gation ; we suppose also that the observations were all made under the 
same conditions, so that they are equally reliable (their relative prob- 
abilities must be determined by the theory). 


II. Case or a Fintre Nomper or OBseERVATIONS. 


Let 2,25 %3,- +» 2% be m observed values of a magnitude whose true 
value is Z, and let z denote an observed value in general. If A is any 
function of z and A,, A,, A;,..., A, are the values of this function for 
21, 25 Zgy-+ +» Zp, respectively, we write 
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A, + A, + Azs+...+ A, = [4], 
as is customary in the theory of errors. Let 2 be the arithmetical mean 
of the observed values, so that 


1 
1 =— [2]. 
(1) a hee [z] 
Let x be the error of an observation in general, x,, 22, 23,..., 2, the 


errors of the given observations, and x, the arithmetical mean of these 
errors, so that 
2=2—Z,2,=2,—Z2,%.=%—Z,%,=%3—Z,...,5=%—Z, 


(2) m => [2] => [2]—Z=%—Z; 


therefore, the arithmetical mean of the errors is the error of the arithmet- 
ical mean of the observations. Let ¢ be the residual of an observation 
in general, — that is, the excess of the observed value over the arithmet- 


ical mean of the given observed values, —and let &, &, &,..., & be 
the residuals of the several observations, so that 
(3) é=2z—%= (4+ Z)— (m+ 7) =zr-NH; 


therefore, the residual of an observed value is the residual of its error ; 
we may call it simply the residual of the observation. Let p, be the 
mean k-th power of the residual, so that 


(4) pm =~ [2s 
in particular, 

p=- (6) =. 01) =1, 
(5) 

pa =2 (== [2] =0. 


The positive square root of the mean square of the residual shall be 
called the mean residual (it is called the mean error in the ordinary 
theory) and shall be represented by pu, so that 

(6) w= Vp, or pap 

The ratio of the residual to the mean residual shall be called the relative 
residual and be represented by s, while the ratio of the mean k-th power 
of the residual to the k-th power of the mean residual shall be called the 
relative mean k-th power of the residual and be represented by p,, so that 
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= fn =F, 
(7) : : 
1 


A little consideration will make it clear that no observations, however 
many, can determine the ¢rue value of the observed magnitude, nor, 
therefore, the errors of the several observations, — they can determine 
only what is equivalent to their arithmetical mean and the residuals of 
the observations. The fundamental problem of the theory of errors 
must, therefore, be modified, and may now be stated to be the determina- 
tion of the probability that the residual of an observation shall lie 
between certain limits according to the law of the long run, that is, in 
general, on the assumption that the distribution of errors over the whole 
range of possible error tends to follow a certain law, as the number of 
observations is increased without limit, and that the proportion of all the 
errors of any (finite) number of observations that fall between given 
limits is the same as if the number of observations was infinite. It is 
evident that the distribution of errors among different intervals in definite 
proportions implies an increase in the breadth of the intervals as the 
number of observations is diminished. The analytical theory of errors 
is based upon the supposition that the number of observations is infinite 
and that the interval for which the proportion of errors is to be deter- 
mined may be taken to be as narrow as we please. 

In connection with this statement of the fundamental problem it may 
be said that the determination of the value that may best be assigned to 
the observed magnitude is not explicitly a part of the fundamental prob- 
lem, although immediately connected with it and perhaps the most 
desirable result of it from the practical point of view, but the criterion 
for the best value is by no means evident. Perhaps it will be generally 
considered that the best value is the probable value, that is, the value that 
is just as likely to be exceeded as not in the long run, so that the prob- 
ability that the observed value shall fall short of the probable value is 4 
and the probability that the observed value shall exceed the probable 
value is 3. Certainly the best value cannot be defined as the most prob- 
able value, that is, the value whose probability is a maximum, because 
it turns out.that the probability is sometimes a maximum for more than 
one value; the probability may even have equal maxima for two or more 
values. 
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Ill. Case or AN Inrinite NumMBer oF OBSERVATIONS — 
ASSUMPTIONS. 


The assumptions upon which we shall base the general theory of errors 
(upon which the ordinary theory has always been based in part, explicitly 
or implicitly) are these : 

I. Possible errors form a practically continuous sequence from a certain 
lower limit to a certain upper limit ; 

II. The probability that the error of an observation lies between x and 
x + dx, where dx is infinitesimal, is (x) ‘dx, where $(x) ts an analyt- 
teal function of x developable by Taylor's theorem throughout the whole 
range of possible error ; 

IIL. The probability that the error lies between given limits is independ- 
ent of the unit of measurement. 

Assumption I forms the foundation for II, and II implies that 
¢ (x + &) is developable according to powers of h by Taylor’s theorem 
for all values of x and x + / that lie within the range of possible error. 
The value of ¢(zx) for any given value of x is called the relative fre- 
quency of the error x. 

The probability that the error of an observation lies between two given 
limits a (lower) and 6 (upper) i . evidently, 


(8) foc -a 


If x is the (algebraically) least possible error and z the greatest possible 
error, we have : 


(9) foe denn 


If we express x in terms of the residual ¢ by (3), representing the lower 
and upper limits of possible residuals by ¢ and é, respectively (which 
correspond to the extreme errors x and z), and writing 

(10) $(z) =o(E+m) =y¥);, 

we have from (9), because 2, is constant, 


t 
_ (11) fv@-a=1, 
£ 
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and the probability that the residual of an observation lies between 
a and 6 is 


b 
(12) f V+ dé. 


The value of y() for any given value ¢ may be called the relative fre- 
quency of the residual ¢, and y (€) « dé is evidently the probability that 
the residual lies between ¢ and ¢ + dé; that is, y (€) * dé is the propor- 
tion of all the residuals that lie between ¢ and ¢ + dé. In calculating 
the mean powers of the residual by (4), the power of each residual is to 
be multiplied by the number of times it occurs, and the sum of the 
products thus formed for all the residuals is to be divided by the number 
of observations ; if we assign all the residuals between ¢ and ¢ + dé to 
the value ¢ itself (or replace them all by their mean), the multiplier 
corresponding to this value is proportional to y (é) * dé, while the whole 
number of observations is proportional to (11). We have then 


(13) Pr = feve * dé, 
é 

and, in particular, by (5) and (6), 

(14) fu@-aan=s, 
£ 

(15) fev@-=n=0, 
é 

(16) fev@-denae 
£ 


Therefore, by (3), (14), and (15), 


(17) furg@yrde= fE+%) VO dan 
z £ 
Because y (é) is essentially positive for any possible residual ¢, formula 
(15) shows that the possible values of ¢ must include both positive and 
negative values, so that §< 0, 0 < , and 0 is a possible residual in 
all cases. 


| 
| 
. 
| 
| 
{ 
| 
| 
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IV. Errect or an InriniresimaL Source or Error. 


To determine the form of the function y (€) we consider the effect of 
the addition of an infinitesimal source of error to the original source 
(however complex that may be), that is, of a source of error that pro- 
duces infinitesimal errors with great frequency and finite errors of even 
moderate magnitude with infinitesimal frequency. We employ the same 
notation for all functions and quantities pertaining to the action of the 
infinitesimal source alone as for those pertaining to the original source 
alone, distinguishing them by an accent, and denote changes due to the 
addition of the infinitesimal source by the symbols for the functions and 
quantities changed with a 6 prefixed. 

A given error z produced by the simultaneous action of both sources 
is the result of a certain error zx’ produced by the infinitesimal source and 
an error x — x’ produced by the original source. The probability that 
the error produced by the combination shall lie between z and x + dx 
is, then, 


[#(e) + 86 (e)]de = f b(e— 2) @) de de, 
v4 


from which follows 
z! 


(18) (2) + 86 (@) = fo e— 2) g(a) det 


For the arithmetical mean of the error produced by the combination 
we have 


z+8e Y 

x, +8" = |x-dx: | o(x—2x')+ d(x’) + de! 

: ie 
z z 

= [xed f $(e—2) $l (e) de! + 2-4@) + B—2-$(0) “be 


z z 


It is to be observed that, in whatever manner these integrals are 
evaluated, x — a must be taken over the whole range of values from x 
to z. The limits of possible error (x and x) cannot be fixed a priori ; 
some may think that large errors or residuals are impossible by the nature 
of accurate observation, but we prefer to regard the limits as determined 
by the frequency function ¢(x), namely, we assume that the value of 
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¢ (x) (which is positive for any possible error 2) is 0 at each limit ;* 
if ¢ (x) is not 0 for any finite value of x, all values of the error are to be 
regarded as possible. The results will show that ¢ (x) = 0 for x = + o,. 
and, in general, ¢ (x) has a very small value for any large value of x. 
We have, then, 


(19) ¢(z)=0 and (zx) =0, 
so that the last two terms of 2) + da vanish, and, on reversing the order- 
of integration, by (17) and (9), 

ry Za 
ty + Bay = f g(a!) del fp (e— 2!) de 
- 


ate! 


=f oe): del- fete) oe) sae 


= { @ + a’) fp! (x!) da! = x) + 2X! 
2 
so that 
(20) Bary = 2 


By (10), as applied to both component sources and their combination, 
we have 


¢! (x') = ¢ (@ + a!) = y' (&), 
¢(x—xe)=O(E—P +m — x) =y (E— f — x’), 
(x) + dp (x) = ($ t+ 34) (E + %) = (YW + BY) (E—%’) = 

v (€ — ao!) + By (E — 20’), 


* The relative frequency y(t) of a possible residual must be positive; if ¢ 
varies continuously from zero, either through positive or through negative values, 
the limit — or — must be reached before y (¢) becomes negative ; but if y (g) passes 
through zero without changing sign as ¢ passes through a certain value, —a value 
that corresponds to a point of contact of the curve represented by y = y (é) with 
the axis of abscissae, —that value of & is not necessarily a limit to the possible 
values; such a value of ¢ could never occur as a residual, although values on both 
sides of it might occur. However, if the conditions under which the observations 
are made change continuously in such a manner that the curve y = yw (¢) assumes a 
form in which it touches the axis of abscissae, this must be the limiting (extreme) 
form, provided a continuation of the change would produce a curve having negative 
ordinates in the neighborhood of the point of contact; namely, on account of the 
permanence of the relation (14), the limits and cannot change discontinuously. 
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for any given value of £, because the arithmetical mean of the error cor- 
responding to the function ¢ + d¢ is 2% + da = 2) + %!. Therefore, 
(18) becomes 


b (E — x) + dp (E — 2’) =fve- & — x!) + y! (é') « de, 
g 
or, when é is replaced by € + 2’, 


(21) y() + EO) =fvE—L)- we): a 
¢ 


for any possible value of ¢ Now, by assumption II, ¢(z) is develop- 
able by Taylor’s theorem and, therefore, so is y (€), that is 


We -> ye He, 
y(é ae iS 1) kl Og ’ 
therefore, by (13), (14), and (15), 


VO+ HO ->. 0 7 f ety (D+ ae! 











->, fue ye % My (é) 


tk 
og* 


am Px MVE), 
— vy () ‘Re 1) k! og* 3 





so that 


(22) By (2) ->, (— ye 


Also, by (13) and (21), 
pt tnnfe-a-[ve—e “y (E) de 
E+E & 


P 
= feat foe) VO) a +E VO BHP VO -% 
é Px la a a8 











STORY.— A NEW GENERAL THEORY OF ERRORS. 175 
Reversing the order of integration, observing that ¢ — ¢’ must take 


all values from ¢ to ¢ and that, by (19), y (2) = 0 and y (é) = 0, we 
have 


retina [Ve ie fer-ve- b) + de 
- 


-fv (e)° w fierervo- dg 
->, c ) fe v (&)° w fe yO ag 
WS, (s)nimeensS,(S)nne 


where (7) denotes the coefficient of z* in the development of (1 + x)” 


by the binomial theorem, that is 


oo ()* orm 
Therefore, a 
(24) Pn “2h ( ) Px Pm for 2 = m. 


In particular, for the smaller values of m, 


ip, = 0, dp,=0 (becausep,=1 and p,=0), 
Sp. = ps', Sps=ps, Sp,= 6 papi! + py 

dp; = 10 ps pro’ + 10 pops! + ps’, 

Sp, = 15 p, po! + 20 ps ps! + 15 ps py! + pe’, ete, 


and, therefore, 


Pa = Spr, py =Sps, ps = dp,— 6p. * Spr, 

ps’ = Sp; — 10 pg * dps — 10 py Sp,, 

Pel = Sp, — 15 py * Spa — 20 ps * Spgs — 15 pa (Sa — 6 pra * Spa) 
= Sp, — 15 py * Spy — 20 pz * Sp, — 15 (py — 6 pq") Spa, ete. 
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Because the right member of (24) is linear in the p’’s and p,,' has the 
coefficient 1 in the expression for dp,, but does not occur in the expres- 
sion of any dp with suffix less than m, p,,’ can be expressed as a linear 
function of the dp’s with suffixes not greater than m whose coefficients 
are polynomials in the p’s. The expression for p,’ in terms of ép,, 
5ps, Sp, and dp, shows that the coefficients of the ép’s in the expres- 
sions of the p’’s are not generally linear. If we substitute the expressions 
of the p’’s as linear functions of the dp’s in (22) we obtain dy (é) as a 
linear function of the dp’s whose coefficients are polynomials in the p’s 
multiplied by derivatives of y (€) with respect to £; therefore, y () ts a 
function of € and the p’s with suffixes as great as 2 alone, say * 

(25) ¥ (€) = F'(E, Pa Par Pas Pos +++) 

This result confirms the statement made above that all theory can 
determine is the probability that the residual shall lie between given 
limits. 

V. Cuance or THE Unit OF MEASUREMENT. 

If we reduce the unit of measurement in the ratio 1 : A, the numerical 
values of all measurements are multiplied by A, so that, for given obser- 
vations, 2, %, &, df, p, (for 2 = &) have to be replaced by Az, Aa, AE, 
dX: dé, A* p,, and, by assumption III, the probability that the residual 
lies between A £é and A é + A + dé with the new unit is the same as that 
it lies between ¢ and é + dé with the old unit ; therefore, 


F(AE, Mpa, A? ps, Mp, .--)*AdE= F'(E, Po, Ps, Pay -)* dé, 
that is 
4 1 
(26) Faé, Mpa, AM ps, Mp, )= xX Fé, Po, Psy Pa, +++) 
for all values of £, pz, 73, Y4,-.- The relations (25) and (26) express © 


that y(é) is a homogeneous function of &, V Pas V/ Pas V Pa, etc. of de- 
gree —1, that is, with the notation of (6) and (7), 


(27) VO= “f ay 


for convenience, we shall sometimes abbreviate y (¢) to y and 


SF (8; Pes Pas pss +++) to f(s) orf. 


* This conclusien implies that every infinitesimal change in y(é), for a given 
value of ¢, can be effected by the addition of an infinitesimal source of error, — 
which ought, perhaps, to be stated as a fourth assumption lying at the base of the 
general theory. 
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VI. DirrerentiaL Equation or f(s). 


Formulae (22) and (24) give the variations of y (€) and p,, due to the 
addition of any infinitesimal source of error ; the mean square p,’ of the 
residual due to such a source is an essentially positive infinitesimal ; in 
determining the mean powers of the residual for such a source by (13) 
we may confine our attention to values of the residual that lie within 
very narrow limits which, in accordance with the remark made in con- 
nection with (15), include 0, so that we need consider only very small 
values of £; the terms of the integral (13) for greater values-of k& are 
infinitesimals of a higher order than the corresponding terms of the 
integral for smaller values of &, so that the integral p,’ for any value of 
k greater than 2 is an infinitesimal of a higher order than p,’. We 
specialize the infinitesimal source of error by assuming p,’=0 for 
3=k* We have then, by (24), 


ép,, = G) Pm-2' Pe for2 =m, 
from which, together with (6), follows 
Pi = Sp. = 2° dp, 
so that ~ Spm = m (m —1)* py_op * Sp, 
while (22) gives 
1 PW) _ Pv.” 
(0 = 4p) SO Oy, 


We denote partial differentiation by the round 9, regarding wy as a 
function of £, pe, ps, Ps». - and f as a function of 8, ps, ps, ps, +++ » 80 
that, by (7) and (27), . 

Rf & OR « 
ofp Os o€?  p® Os? 

The variations denoted by 5 were all taken for a definite value of £; 
therefore, by (7), (27), and the formulae just written, 


é é bu 
=$§(-)=—-—°d=-—s-:—» 
ds (<) 5 Op 8 


d oF 8 
y= Tn ea Se 





an 





* This assumption does not in any way affect the generality of the results 
obtained, because p,’ enters into y (¢) only through the variations of the unac- 
cented p’s. 

VOL. xL.— 12 
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2 8 
Y= suv) =n + y-m= [Sts] *, 
dp, = k (k — 1) * ppsp*: 8p, 
ne Px _ Px ene mie 2 ee a op, 
t= 8 (2) = Br — eB bn = (EE —1) mak] 2 


therefore, 


aE tet 2 tnt bp + bos + 


=H + Dag, 
={-.% +2, [- (k—1) prs — hor] |. ‘ 


Equating the two expressions for $f thus found, omitting the common 


factor Su » we obtain 
B» 


(sy = SF 4 t=, [FE —D)me—kod Zo, 


3s a + 
a differential equation for f which may be written 
(29) 2s 5 (a+ sf\}=2: Tf, 
where 
d 
(30) T= [EE ma— tod 5 


is a differential operator that acts upon f only through the p’s with suffixes 
as great as 3; the first part of that term of 7’ for which k = 3 vanishes 
because p, = 0. We have omitted the upper limit to the sum in 7’ be- 
cause that depends upon the number of p’s involved in f, a number as 
yet unknown. 


VII. Tae Operator 7. Tue P’s anv S's. 


The operator 7’ is so important for the theory that it deserves par- 
ticular study. By a constant we shall henceforth mean a quantity that 
is independent of the p’s with suffixes as great as 3 as well as of s. If 
u and v are any two operands and C is any constant, we have 
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T(ut+v)= Tut 7'r, 
(81) T:Cu=C:Tu, 

T-uv=u‘'Tv+ov° Tu. 
If u is a function of the p’s with suffixes not greater than m, 7'u is also 
a function of the p’s with suffixes not greater than m; if u is a poly- 
nomial in the p’s, Z'u is a polynomial in the p’s of the same degree as wu; 
in particular, if u is a linear function of the p’s with suffixes not greater 
than m, 7'u is a linear function of the p’s with suffixes not greater than 
m Its em, 
(32) 7’ pm = 4m (m — 1) * pms — 3 ™ * Pw 

Therefore, if u is a linear function of p,, and the p’s with suffixes less 


than m, where 3 = m, 7" u is a linear function of p,, and the p’s with 
suffixes less than m in which the coefficient of p,, is the coefficient of p,, 


in « multiplied by | — ;) , where ” is any positive integer or 0; more 
generally, if G(Z’) denotes any polynomial in 7’ (linear function of 
powers of 7' with constant coefficiénts), G(7')* wu is a linear function 


of the p’s with suffixes not greater than m in which the coefficient of 
pm is the coefficient of p, in u multiplied by G (- 3)3 in particular, 
G(T) : u will not involve p,, if G(7’) contains 7’ + 5% a factor, and 
only then. Formula (32) holds also for m = 0, 1, and 2, giving 

Tp =9* pp =90, Ty =—tn=9%, Ter =m —p=9, 
because po = pz = 1 and p, = 0, by (7) ; therefore, 

(T+ %) po = 9, (7+ 4) mr =0, (7+ 9) r= m= 15 
so that G(7')* p» does not involve p,, if G(7') contains the factor 
T+ 3 , where m is 0 or any positive integer excepting 2, and G(7’) « p, 


does not involve pz if G(7’) contains 7 as a factor. Of course, poly- 
nomials in 7’ are commutative with each other, so that it is of no conse- 
quence what position in a product G (7’) the annihilating factor occupies. 

It is convenient to employ the usual symbol y” to denote the product 
of n factors diminishing by successive units from y, whatever meaning y 
may have (it may be a symbol of quantity or of operation), so that 


33) yy =y¥ (y—1) Y—- 2) Y—-3).-- Y—2+)) 


ee 
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if n is any positive integer; in consequence of this definition, 





%=y, 
(34) yay (y—n), 
yo) 
f= , 
ee 
for 1 ='n; we take the last formula as the basis of the definition of y”’ 
for n = 0, — 1, — 2, — 3, . . ., successively, and thus obtain 


ad 1 
wah £"=GEHNGFDYt9 Gn) 


for any positive integer m (that is, for any negative integer —n). With 
this notation we have 


fy = : (y ie Yate ae (— 1)" ° (— ytn— 1), 3 
(y +n— 1)” —_ (— 1)" ‘ (— y)™, 





(35) 





(36) jn™=n!, nv? =——_, 


(7) n'? re a r) (; _ i) 

. => ° “++ ° ’ 

\\e a! t ti—1 

where each of the letters m and 7 represents any positive integer or 0; 
also, 





((2%)!= 2*-a!@— 4)” 
(37) \ . — 9+. SIs +1 
(2¢+1)!=2 ale + 4yer, 
because, evidently, 2‘7! is the product of the even numbers from 2 to 2 /, 
2‘(¢ — 4) is the product of the odd numbers from 1 to 27 — 1, and 
2'+1 (@¢ + 3) is the product of the odd numbers from 1 to 27+ 1. 
Furthermore, we have 


(y — 1)" = (y— 1)! = (yn) = y= (y— 
= — ny) + = yD — 2 - yr, 


(38) : 
=> Hr ney, 
0 


where n is any positive integer or 0. 
From (32) it is evident that any polynomial in 7’ changes p with any 
even suffix into a linear function of the p’s with even suffixes and p with 
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any odd suffix into a linear function of the p’s with odd suffixes. There- 
fore, if «is a linear function of ps,, and the p’s with even suffixes less. 
than 2m, where 2 = m, (7'+ m) u does not involve p.,,, (7+ m)” u 
involves neither pg, DOT poms, (7'+ m)® wu involves neither pon, Pen—2s 
NOF Pom—4, and so on, and (7'+ m)— u is a constant (or linear function 
of po and pe); therefore, by (36) 


ss ‘hate m)™*) y = (T+ m)™ (T+ 1) 
09) =(T+1)7T:(T+m"u=0; 


that is, (7’+ m)™*” annihilates every linear function of the p’s with: 
even suffixes not greater than 2m, where m is any positive integer or 0» 
(this is true for m = 0 or 1 because the operator contains 7’ as a factor). 
But (7'+ m)*» annihilates no linear function that involves a p with: 
an even suffix greater than 2m or an odd suffix as great as 3; namely, 
if p, is the p with greatest even suffix or greatest odd suffix (other: 
than 1) that is involved in the linear function u (with a coefficient dif-- 


ferent from 0), the coefficient of p, in (7'+ m)"* wu is the coefficient 
(m+1) 

of p, in « multiplied by (m — 5) » whose factors are all negative if 
kis an even number greater than 2 m and all odd multiples of 4 if & 
is odd. 

Similarly, if « is a linear function of p2,,,; and the p’s with odd suffixes 
less than 2m + 1, where 0 = m, 
(40) (7+ m + 3)" u = 0, 


so that (7’+ m+ 4)*” annihilates any linear function of the p’s with 
odd suffixes not greater than 2m + 1; but (7’+ m + 4)*" annihilates 
no linear function that involves a p with an odd suffix greater than 
2m +1 or an even suffix. 

Applying (39) and (40) to po, and pon4;, respectively, we have 


(T+ m)™* pom = 0, (T+ m + 4)" pons = 9, 
so that, by (34), 
T (T+ m—1)™ pom = — m* (T+ m— 1) pons 


T(T + m— $)™ poner = — (m + $) * (T+ m — $)™ pants 
or 


(41) T Qen = — M* Qoms T Qom1 =— (m+ 4) 5 Qam+ts 
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if we write 
(42) (2+ m—1)™ pan = Qem, (T+ m — $)™ panes = Qomtt j 
and, evidently, formulae (41) hold for any constant multiples of Q,,, and 


Qensi» Tespectively. 
Now (32) shows that 


(7+ m) * Pam = 2m (m — $) * pans 
and, by continued application of the same formula, 
(48) (T+ my” pay = 2 m! (m — 3) * pan 
for 0 = k = m—1;; in particular, 


(T+ m)'""—") ps, = 2-1 pl) (m oe 4)" ps 
and 
(44) (T+ m)™ poy = 2" + m!(m— 4)™ pp 


because, by (34) and (36), 


m!=m™ =m) +1, 


(m — $)™ = (m — 9)" +f, 
(7+ 1) ps = po = po; 
so that (43) holds also fork = m. Then, by (42), (38), (48), and (36), 


Q =D Dim (7+ my pm 


= 2" (m — 4) > m >, 2 (7 ) FGF 
for 0 = m. 


Similarly, we have, from (32), 
(7+ m + 4) * pomsr = 2m (m + $) * pons 
and, by continued application, 
(48) (T+ m+ 4)" pangs = 2m (m+ 8) * pan ans 
for0=k=m. Therefore, by (42), (38), (46), and (36), 


(45) 


Qanis -2, (— 1) m9 (T+ m+ 4)" pangs 


(47) 
= 2" (m + 4)" + m >, (— (7 ) gar et pai GT pa yp 


for 0 2 m. 
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aoe ree De 


we have, by (45) and (47), 
Qon = 2" (m — 4) m! P,, and 
{ omy == Lt? (m + $) "9 m! Paar, 
“and, because P,, and P,,, are constant multiples of Q,,, and Qon4:1, 
by (41), 
(50) T Px», =—m* Py, PD’ Pmgs = — (mM + $) * Ponts 
for all values of m. For the P’s with the smaller suffixes we have, from 
(48) and (7 )» 
(51) PB Pp=P,=0, Pr=—te, P=—l+tn, 
Po=—fpstrses, Po=—2+pi—rspe ete; 
so that equations (50) give 
(52) TP,= 0: Py = 0, TP,=-—}P,=0, TP,=—P,=0. 
Because P;,, is a linear function of the p’s with even suffixes up to 2m 
in which the coefficient of p,, is different from 0, and P,,,,, is a linear 
function of the p’s with odd suffixes up to 2m + 1 in which the coefficient 
Of pon41 i8 different from 0, p,,, can be expressed as a linear function of 


the P’s with even suffixes up to 2m and p2,,1 a8 a linear function of the 
P’s with odd suffixes up to 2m+ 1. In fact, if we multiply P,,, or Pons 


by (— 1)” (.) and sum for all values of m from 0 to k, we obtain 


S.cm(f)a “Se eee oa )(7) 
is Fe- 


>: 1 (7) Pat ->; CU ea sree ke ony?) 


If we write 
; 


(48) 





(49) 





= pei Es + pe 
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(cro) 
-—9()Rcr('s) 


=Oifi<k and 
= (— 1) if +=&k; 
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because 


Ser) 


therefore 


Pa = 2* (k — p> ce 1)" a) Py, for 0 2 k, 
0 


(58) ; . 


pas = (E+ NS Ky ( 1) Panes for OZ 
m 
\ 





0 
As a consequence of (50) we have 


(54) G(1)P,= @(-3)-P,, 


if- G(T) is any polynomial in 7. The ease with which the result of 
applying any polynomial in 7’ to the /’s can be obtained by means of this 
formula makes the introduction of the P’s in place of the p’s very advan- 
tageous, We shall henceforth assume that all functions of the p’s have 
been expressed as functions of the P’s by (53), so that the degree of any 
polynomial in the P’s with suffixes as great as 3 is the same as that of the 
corresponding polynomial in the p’s. Equations (50) show also that the 
operator 7; defined by (30) in terms of the p’s, can be expressed in terms 
of the P’s thus: 


9 
(55) fax > Pap 
3 


We assign to each P a weight equal to its suffix (leaving out of account 
P, and P,, which are 0 and, therefore, nowhere occur) and define the 
weight of any product of P’s to be the sum of the weights of its factors ; 
any constant may be regarded as a multiple of P) = 1, and, therefore, as 
of the weight 0. The weight of any such product as Ps" P,% P,% .. . is 


(56) ‘ 8n,+ 4n,+5n,+...= 0, 
say, so that, by (55) and the last formula of (31), 
(57) TP P™P®™...=— 3 PP Po. ., 
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and, evidently 
(56). O(T)* RPP MP. ..= e(- 5) PPP rR ess 


where G(7') is any polynomial in 7. More generally, if 1 is any 
isobaric polynomial of weight w in the P’s, that is, any polynomial whose 
terms are all of one weight w, 


(59) G(T): = o(-$) =. 
Therefore, if the coefficients of II are not all 0, G(7’) - II is or is not 0 


according as G(7') does or does not contain 7’ +5 as a factor. In 


. particular, (Z7'+ m)™*” annihilates every polynomial in the P’s whose 
terms are all of even weights not greater than 2m and annihilates no 
polynomial of other description, while (7’+ m + 4)™*” annihilates every 
polynomial whose terms are all of odd weights not greater than 2m + 1 
and no others. By means of formula (59) any polynomial in 7’ can be 
applied to any polynomial in the P’s; namely, any polynomial in the 
P’s can be expressed as the sum of isobaric polynomials of different 
weights and the polynomial in 7’ can be applied to each isobaric part 
by (59). 

The P’s can be determined from the function f(s); namely, from (7), 
(18), and (27) we have 


Z ; 
k 
(60) p= B= fio a= st f(s): ds for0 Zk, 
§ 


8 


and obtain, by the substitution of these values of the p’s in (48), 


(61) Pa= fs “F(8) ds, Pamss a Se “f(8) « ds, 


2-300 (Tega 


eux ->- A ces 


where 


(62) : 





for 0 = m. For the S’s with the smaller suffixes we have 
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SH=1, S=s, 2 = 1—s’, S;=s— 38, 
S=1-28+ 48, S=s— gs + p85, 
S, = 1 — 38? + s* — y5 58°, ete. 


VIII. Tae Form or f(s). Tue R’s. 
To simplify the integration of equation (29), we put 
f(s) =0(8) 4, 
where wu is an integral of the equation 


Ou 
| eines u= 0, 


say 
8? 
wae 3%’ 
where e represents the base of natural logarithms, so that 


(64) f(s) = (8) e-3. 
Writing @ as an abbreviation for 6 (s), we have 
of [00 .. it, 
os <a i ale : 
OF 30 


=| ge 2" cea: “3 
so that (29) becomes 


E& 3 ol o-$ a 0-T6-e-8? 


because e-3 is independent of the p’s (and P’s) ; therefore, 


(65) ae a= 2 T6, 
which is a differential equation for 6 (s).’ 

The assumption II that ¢(x) is developable by Taylor’s theorem 
‘throughout the range of possible values of x implies that y (¢), f(s), 
and 6 (s), are similarly developable; we may therefore put 


(66) 6 (s) => 7 si, 
0 


where A, (for each value of 7) is a function of the P’s alone. This de- 
velopment implies that we shall approximate more and more nearly to 
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the true form of 6 (s), even if its complete expression is an infinite series, 
the more terms of (66) we take into account. Taking any finite num- 
ber of terms, say for values of ¢ up to «x, substituting in (64) and the 
resulting expression for f(s) in (61), we obtain 


(67) P.. ->, a fe R,, « ds, 
0 s 


where 
(68) R,, = S,*e7% for 0 = m. 


We see from (67) that, to the degree of approximation that corresponds 
to x + 1 terms of (66), the P’s are linear functions of the x + 1 A’s, 
whose coefficients are definite functions of the limits s and s alone, because 
the S’s and, therefore, the R’s are independent of the P’s. Therefore, 
not more than «x + 1 of the P’s (including P,, P,, and P,) are linearly 
independent ; the solution of such equations (67) as correspond to linearly 
independent P’s gives the x + 1 A’s as linear functions of these P’s, 
whose coefficients are definite functions of s and s derived from the inte- 
grals involved in (67); these expressions of the A’s will involve a certain 
number of arbitrary constants if the number of linearly independent P’s 
is less than x + 1, but that is of no consequence. If the P’s are sub- 
jected to infinitesimal variations, —say SP,, is the variation of P,,,— 
the A’s, f(s), and @(s) will undergo variations, which may be repre- 
sented by 84, (for A,), 8f(s), and 50(s), respectively. Then, by (66) 
and (64), 



























(0) => “4 TT" df (s) = 86 (s) -e-? 


and, by (61) and (68), because f(s) = 0 for each of the limits s and s, 


SP, -fs. f (8) ° un fr 80 (s) - 0S fer, ds; 


that is, the relations between the 5P’s and the 5A’s are the same as if 
the coefficients of the A’s in (67) were constants (whether s and s are 
functions of the P’s or not); therefore, in determining the A’s as func- 


tions of the P’s that satisfy the differential equation (65) we may regard 
them as linear functions of the P’s with constant coefficients. The 


necessary and sufficient conditions for these linear functions are those 
implied in (65) together with 
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(69) A,=0 forx <i 


(to the degree of approximation to which we restrict our calculations). 
Substituting the expression (66) for 6 (s) in (65) we obtain 


i ae ee at 
ars (F— 2) ay (@— 1)" = > a7” 


» 3 [ 4a ihe 4, | 5 oe 


0 
i 
for every possible value of s, where the first term of the multiplier of # 


that is 


is to be omitted for 7 = x and 1 = x — 1 and the second term goes out 
of itself for i= 0. Therefore, the necessary and sufficient conditions 
for the satisfaction of (65) by (66) are 


A= 2 (745) °4, for0Zi=n—2 


and 

(7+5) *A,=0 for t=« or ex—1; 
that is é 
(70) 4,=2(7+°S*)- 4, for? Zi Ze 
and 


(71) (T+5) 4=0 


for A, with the greatest even suffix and A, with the greatest odd suffix 
that occur in (66). On replacing ¢ in (70) by 2¢ and 24 + 1, respec- 
tively, we have, for the A’s with even and odd suffixes separately, 


Ay =2(T+¢—1)* Aggy for 1Ziz5, 


Aga, = 2(T+i—})'A,, for 1S 
and, by repeated application of these formulae, 
(72) Ay, = 2(T+i—1)+A, for 1Z¢= 
and 
(783) Any = 2(P+i—}): A, for 1S¢S ne 
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that is, the A’s with even suffixes are all derived from A, and all the A’s 
with odd suffixes from A, by the application of definite polynomials in 7. 
We may, then, determine the forms of the A’s with even suffixes and of 
those with odd suffixes separately ; that is, we may suppose our calcula- 
tions restricted to a degree of approximation that corresponds to a devel- 
opment of 6(s) in which occur even powers of s up to s?* and odd powers 
up to s?8+1, without assuming any particular relation between a and f. 
Then (72) holds for 1 = t S a, (78) holds for 1 = ¢ = B, while (71) is 
expressed by 
(T + a) * Ao = 0, (7+ B+ 4) ° Ags = 0, 

that is, when Az, and Aogi; have been determined by (72) and (73), 
respectively, 

(74) (T + a)@t) + A, = 0, (7+ B+ 44: 4,=0. 

The only conditions imposed upon A, and A, so far are, then, that they 
shall be linear functions of the P’s with constant coefficients that satisfy 
the equations (74). The only dinear functions of the P’s that are of 
a given weight w are constant multiples of P,,, while, by (54), 


(a+1) 
(T+ ayer Py = (a3) * Per 


(T+ B+ yer Py = (p—2S")""-p,, 


therefore, in order that A, (as a linear function of the P’s) shall satisfy 
the first equation (74) it is necessary and sufficient that the suffix w of 
any P involved in it shall satisfy the equation 


(a+1) 
(--1)""=4 


that is that w shall be an even integer not greater than 2a; and, in order 
that A, (as a linear function of the P’s) shall satisfy the second equation 
(74) it is necessary and sufficient that the suffix w of any P involved in 


' it shall satisfy the equation 


w — 1\(6+1) 
Gar lien 


that is that w shall be an odd integer not greater than 28+1. We 
may, then, write 


a 8B 
(75) Ay =>, Cy Px, A, =>, % Pris 
0 0 


where the ©’s are constants. 














ees. AeA. aC 


Re eee mia 
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On applying (72) and (73) to the first and second of equations (75), 
respectively, we obtain, by (54) and (36), 


A=), Ce Py 2 (i—k— =>, Oy Py, (—1)' 258, 
Om. “, eae? 
ee =>. Cnnss Panyr* 2(¢—k—1) =>, Onna Panyr* (— 1 25, 
0 0 
Substituting these values of A,, and A,,,, in (66) we have 


g BO si 
6 (*) ->, Cn Pa ra Bie aa 
+, Corrs Pass > yon Yaz + api 


or, replacing (27)! and (2% + 1)! by their values in (37), observing the 
last line of (36), and employing the notation of (62), 


a B 
(77) 6 (8) = > Ci Pu Su + > Cres Paver Soesrs 
> 2 >» ak+1 4 ae +1 
again, by (64) and (68), 


a B 
(78) F (8) =>,.% Py, Re +>, Cons Props Rowyr- 
0 0 


Substituting this — of f(s) in (61) we obtain 


=>, Cu Px J sate aS, Ore41 Pores Jstawe , 
> Coz Pox" fate uw, Orns Pox fata ds. 


(79) 


IX. Tae Coerricients or f(s). Conprtions ror THe Limits. 


In all that contributed to the determination of the form of f(s) as 
given by (78) no assumption has been made with regard to the values of 
the P’s with even suffixes not greater than 2a and the P’s with odd 
suffixes not greater than 28 + 1 [excepting that, by (51), P, = 1, 

' 
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P, = P, = 0], —in fact, they have been assumed to be entirely arbi- 
trary; therefore, the first of equations (79) for 0 =m <= a and the 
second of these equations for 0 = m = f must be satisfied independently 
of the values of the P’s involved; that is, for any one of the values of m 
in question, the coefficient of aay P in the right member of the corre- 
sponding equation must be 0 excepting the coefficient of the P that 
occurs in the left member, which must be 1. Therefore, 


5 
Con {Sim Ban * dt = 1 for 0=mZa, excepting m=1; 
8 


(oy) 2. 
Cants ff Sinss Rania * a8 = 1 for leamZB; 
g 


J Sun Bay de = 0 for OF m= a,0ZkSakSm; 


ry 
J SowBuys * do = 0 for OSmea0EkeBPB; 
(81) /s 


8 
% 


3 
J Sines Bass * ds = 0 for 0Sm=R,0ZkKSRESm. 
A 


From the first, second, third, and fourth sets of equations (81) we ought 
to exclude cases for which & = 1, 0, 1, and 0, respectively, because 
P, = P, = 0 and the terms corresponding to these cases do not occur in 
(79) ; but, by (68), 


5 = s 
(82) JSoBivde= [8.8 oF de = [SR ds 
& & & 


so that the cases that are to be excluded from any one set of equations 
(81) will occur (with an interchange of m and &) in other sets, excepting 
only the case of the integral 
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3 ¥ 
fSikiods= [sas 
& s 


which does not occur as the left member of any equation. The four sets 
of equations (81) can be represented, by virtue of (82), as the one set 


(81’) [Sma = 0 
& 


for any two different numbers m and & from among the even numbers 
not greater than 2a (including 0) and the odd numbers not greater than 
2B + 1, excepting m= 1, k= 2, and m=2,k=1. The equations 
(80) determine all the constant coefficients CO that actually occur in (78) 
(C, and C, do not occur), while equations (81) express relations that 
must exist between the limits s and s for any observations that conform 
to the theory. 

To determine the values of the integrals that occur in (80) and (81), 
we consider the functions S and R defined by (62) and (68). We have 


= (0) GP 
=», (FI saa 
=-2 >, ("7") saqg— pao 
> Xe I rca 


that is, 
(83) oA = — 2m Sa, Ht = 8, 


for0 =m. We have also, by (36), 


f{m—1 g7i+2 
= Soma eer, 7 eer 1) ( i \ = (i + 4)e0 
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(84) = Sana — 8 * Sani for 0 = m, 


m i+1 
S* Som =>- 1)' (*) [2m +1— 2(m — yey opre 
0 
m—1 


: pe | 2+1 
= (2m+1)° Samy — 2 * an C »(" PF \= Gz 4)» 
= (2m+ 1) 4 Som+1 —2m: Som—1 5 





that is, 


2 
(85) Sanit = 5— 


1 
2m+1_ 2m +1 
From (68), (84), and (85) we obtain 
Rom = Rena — 8 * Roms for 0 = m, 
(86) 


Som—t1 + ————_ *$° Sen for 0 = < m. 


2m 1 rors 
Rimi =a Bm toe +1 ‘s’' R,, for 0 =m, 


while from (68), (83), (84), and (85) follow 


Rom ae 
ee 


(87) =—(2m+1)* Bangi; 


s? s? 
onan = (Som — 8 * Samgi) @ ; Soe Gs F 


— (2m * Sony + 8° Syn) O72 = — (2m + 1) Seng O77 


“ — Renae . 
From (83), (87), and (88) follow, after simple reductions by (84)-(86), 


S fm * Sams Bar — + Son Bans} = (0B) * Sax Be 
é {Soms1 Rox — Som Roxyi} = 2 (m — k) * Somer Rory, 
5 42m Sams Banas + Som Rus} = (2m — 2b —1) * Soa Rens 


Therefore, letting cr denote the excess of the value of the enclosed 


function of s for s = 8 over its value for s = s, we have 


(89) [Sia Ra ds = =" (Sas Ral, — =~ [Sm Renal, 


VOL. xt. — 138 
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for ks > mM, 
() 


fen Ba: =e 


5 (Samis Buel, — gq Sm Reel 


for k Sm, and 

(1) 

fata ds = ss [Sins Passs], + ee [Sim Pax], 
Also, from (83), (87), and (88) follow 


S (Sin Roms) = Sam Ram — 2m * Som 1 Romy 


© (Sass Ran) = — (2m +1) * Sines Ramis + Sim Ran 


so that 


(92) fue ds = [Sin Raya), + 2m * fsoeter ds, 


(93) it Oe Seay Same Re}, + 55 rif SR ds. 


By alternate applications of (92) and (93) any lnsegeal of the type 


& 
fs. R,, * ds can be reduced to a linear function of expressions of the type 
8 


[S, Ral for values of & descending from m to 1 and of the integral 


(94) [aman fota 


X. Case IN WHICH ALL VALUES OF § ARE POSSIBLE. 


For s = — @ and s = + « (that is, if all values of the residual are 
possible), the integrals in question have simple numerical expressions. 
Namely, whatever integer may be, 


2 
se 7=0 fors=+; 
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therefore, 

(95) R, = 9 fors=+0; 

therefore, equations (81) are all satisfied, by (89)-(91). Also, (92) and 
(93) become, for this case, 


to +@ 
[Sm Rom * ds = 2m f Sim Ron * ds, 


+o ey 
f Somes Ramis * dt = oz f Sim Ban * ds 


by alternate applications of these two — we obtain 


boo 
f Son R,,, * ds = af Sam—2 Ramo * a8, 


boo 
J Sanss Bamis * do = fete ds; 


and again, by repetitions of these, 


+e 
JSn Punto omfn ds, 


+20 
. m! 
[ a ee ee 5, Bde 


Now, by (93), (63), and (68), 
re ; as s? 
for-aafgn an e *+ds= V2rn, 


where w represents the ratio of the circumference of a circle to its diame- 


ter; therefore, 
+o 
! 
(96) JS Ban te = VER, 


—o 


— ; 
(97) f Ssnti Romyi * 4 =m wep ¥ v2n, 








| 
| 
| 
: 
Hy 
| 
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so that, by (80), 


— 4 
(38) G= ci © # for 0 = m = a (excepting m = 1) 





V274 m 
and 
ea 3. 
(99) Cons ie Re gy for 1 <me B. 





It will be observed that C,,, is less and C,,,,; is greater than oe 
27 


for very large values of m we have, by Stirling’s formula, 
cw. See aN Mem 
Ll / m T 
so that C,,, diminishes and C,,,,, increases without limit as m increases 
without limit; this fact will, perhaps, to a certain extent offset the con- 
clusion that might be drawn from (61) that the values of the P’s with 
odd suffixes are likely to be less than those of the P’s with nearly the 
same even suffixes (on account of the presence of the odd functions S 
under the integral sign in the former and of the even functions S in the 
latter) in estimating the relative importance of the terms of f(s) as given 
by (78) that contain R’s with odd suffixes and R’s with even suffixes. 
For the C’s with the smaller suffixes we have 


1 3 





(100) Con 


qg=— = Ss » = —  _=» ete. 
(101) ae || eee pa .. 
3 15 35 





0; = —— 9 C. = ———_» CO, = ————_; te 
. aes... O00e \" tie 
so that, by (27), (78), (68), (51), (63), and (101), 


0 = 1/0 = = [1-30 —$p) (1— 2a + Bh) 
—$§(1—dpat sooo) (1 — 358? + st — yys®) + ete. 


— 4p, (8— 48) — 4 (ps — Yo ps) (8 — 38° + ys 58) + ete.Jo-F. 





(102) 


XI. AppLicaTION TO ANY ParticuLtaR Case. 


It seems reasonable to suppose that, in any case (whatever the values 
of s and 8) f(s) will have small values for values of s beyond the limits 
of possible values, so that the extension of the limits to — 0 and + x 
will not make much difference in the values of the integrals by which 
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the C’s are determined; we may, then, take the values given by (98) 
and (99) as approximate in all cases. 

The restriction of the development (78) of f(s) to a finite number of 
terms is justified only by the fact that the P’s that would occur in later 
terms are all negligible, that is that 
shag for a<m, 

Pom =9 for B<m, 

to the degree of approximation to which the calculations are restricted. 
In the special case for which s = —o and s= + o these equations 
(103) are exactly satisfied by (89)- (91), because the suffixes of the P’s 
in (103) are different from those of any of the P’s that occur in the 
right members of (79). 

In applying the general method to any particular set of observations 
we have, then, to calculate, by (1)—(7) and (48) (from the observations), 
the values of z, the residuals (€ = z — 29), the p’s, the p’s, and the P’s 
with so many successive suffixes as to satisfy ourselves that the P’s with 
suffixes greater than a certain number are all negligible; the greatest 
even suffix and the greatest odd suffix of the P’s that are not negligible 
are the values of 2a and 2 8 + 1, respectively ; the values of the P’s 
that are not negligible and the approximate values of the corresponding 
O’s given by (98) and (99) being substituted in (77) determine an 
approximate form of 6 (s); the positive and negative roots of 0(s) = 0 
that have the smallest absolute values are approximate values of s and 3, 
with which more accurate values of the Q’s can be calculated by (80) ; 
substituting these values of the @’s in (77) we obtain a more exact form 
of 6 (s), from 6 (s) = 0 more exact values of s and s, and from (80) more 
exact values of the C’s; and so we continue until the variations of the 
values of the C’s are reduced to negligible quantities; the final values 
of the C’s being substituted in (78) give the most exact form of f(s) 
and, by (27), the most exact form of y (¢) that can be derived from the 
observations. The closeness with which the conditions (81) are satisfied 
by the final values of s and s may be regarded as a test of the goodness 
of the observations, that is of the closeness with which the residuals 
follow the law that they follow most closely. 

It may be observed that there is no case for which a = 1; more 
properly, the cases for which a = 1 coincide with the cases for which 
a=0. 

If Panui = 0 for 1 = m, f(s) as given by (78) is an even function of 
s and the relative frequencies of positive and negative residuals that have 


(108) 
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the same absolute magnitudes are equal, that is, positive and negative ‘ 


residuals are equally probable. 
The calculation of the true values of the C’s by (80). (92), and (93) 


will be facilitated by tables of the values of S,, 2,,_, and of the integral 
(94). The values of S,, #,,_; will be required for such values of m as 
are suffixes of the P’s that are not negligible, which will be small, un- 
doubtedly, in all actual cases ; tables of S,, R,, for m = 3, 4, 5, 6 will 
probably be sufficient. For the values of the integral (94) the ordinary 
tables of 


(104) &(y) = “ f e-? + dt 
0 


can be used ; namely, 


4 2 
fetasva oa = 77-0 
0 0 
so that, observing that ® (— 7) = — ® (7), we have 


oo frtanr [o(5)-4(5)] 
Vi [(5)*@)} 


XII. Propasre VAvves. 

We have defined the probable value of the observed quantity to be that 
value that is just as likely to be exceeded as not ; if the probable residual 
é, and the probable relative residual s, be taken to correspond to this 
probable value z,, ¢, and s, will be determined by the equation 


é 8, 
(106) [v@-a=[ro-ans 
3 £ 


and z, by 

(107) % = % + &, 

in.accordance with (3), so that the arithmetical mean is (in general) not 
the probable value even of a quantity directly observed (as it is assumed 
to be in the method of Gauss). 


where s and — s are both positive. 
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XIII. Tue Orpimary Case. 
In the ordinary theory it is assumed that 
(108) P,=90 for 324, 
so that the only P regarded as not negligible is P, = 1, which being 
substituted in (77) with the approximate value C, = sig by (101), 


V 20 
gives, by (51) and (63), 
1 
G(s) = ORS = —) 
( ) 0 Py So oe. an ’ 
6 (s) = 0 has, then, no real finite root, so that 
s=—, s=+o0 
and (80) gives again the same value of Q, as (101), which is, therefore, 
the most exact value possible on the assumption made. We have fur- 


ther, by (78) and (27), 
(109) f(s) = 
while (60) gives 


1 


48 a 4 
—e y= —e 2; 


V2 x 


oo 
=e s&§-e %+ds 
gs week . 


that is 

(110) Pom = 9 for OSm 

and 

(111) pam = 2" (m@—F)™ OF Pay = 2” (m — BY™ + em 

for 0 = m, so that, far from being negligible for the greater values of 
M, pom increases beyond all limit as m increases. 


XIV. Two Orner ParticuLar Cases. 
For the purpose of illustrating the variety of forms the general theory 
may take in particular cases, and especially its possible deviation from 


the ordinary theory, we have plotted the curve represented by y = : F(s) 
B 


(the probability-curve) for a = 0, 8 = 1 and for a = 2, 8 = 0, using the 
approximate values of Q), O;, and C, given by (101), for such different 
values of ps and p, (in the two cases, respectively) as to show all the 
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essentially different forms for these pairs of values of a and 8. Figure 1 
shows the locus of the equation 


1 _f1-ga(s—gje™ 
uV2n ee 


for pg = 0, 1, 1.648, 2.2, 3, 5, and Figure 2 the locus of the equation 


(112) y= 


1 oe 
(113) Pa 8 at ee eee -d 


for py = 0, 1, 1.4, 2, 3, 4, 6, 7, 11. The centre of each figure corre- 
sponds to s = 0 and the two extremes to s = + 4; the horizontal line 
is the axis of abscisse. It will be seen that the ordinates of all the 
curves are practically 0 for values of s whose absolute magnitudes 
exceed 4. The curve for which ps = 0 in Figure 1 and that for which 
p, = 8 in Figure 2 coincide with the ordinary probability-curve. The 
curve for which p; = 3 in Figure 1 and that for which p, = 7 in Figure 2 
touch the axis of abscissae ; these values are, by the footnote on page 173, 
the extreme values of p; and p,, respectively; each figure contains one 
curve that corresponds to a value beyond these limits, which shows the 
discontinuity of s or s a3 the curve passes through the form that has 
contact with the axis. While p, is essentially positive, as a mean even 
power of the relative residual, p, may be negative; if the sign of the 
value of ps; that corresponds to either curve in Figure 1 is changed, the 
curve will be reversed, right and left. In the cases represented in 
Figure 1, the maximum relative frequency and the probable value of the 
observed quantity correspond to negative values of s; in the cases rep- 
resented in Figure 2, the maximum relative ffequency corresponds to 
s = 0 for some values of p,, while there are two maxima for the smaller 
values of p,, the probable value of the observed quantity corresponding 
always to s = 0. 

Because s = 1 corresponds to ¢ = uw, that is to ¢* = p., where p, is 
the mean square of the residual, it is evident that at least one of the 
limits s and s must be considerably greater than 1 in absolute magnitude, 
in any case; but the curves plotted show that very little of accuracy is 
gained by considering values of s beyond + 4 in the cases represented 
by them. 


Crark University, Worcester, Mass. 
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